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ABSTRACT 


Extant chemical evolution models underestimate the Galactic production of Sr, Y 
and Zr as well as the Solar System abundances of s-only isotopes with 90<A<130. 
To solve this problem, an additional (unknown) process has been invoked, the so- 
called LEPP (Eight Element Primary Process). In this paper we investigate possible 
alternative solutions. Basing on Pull Network Stellar evolutionary calculations, we 
investigate the effects on the Solar System s-only distribution induced by the inclusion 
of some commonly ignored physical processes (e.g. rotation) or by the variation of the 
treatment of convective overshoot, mass-loss and the efficiency of nuclear processes. 
Our main findings are: 1) at the epoch of the formation of the Solar System, our 
reference model produces super-solar abundances for the whole s-only distribution, 
even in the range 90<A<130; 2) within errors, the s-only distribution relative to ^^°Sm 
is flat; 3) the s-process contribution of the less massive AGB stars (M<1.5 M©) as well 
as of the more massive ones (M>4.0 Mq) are negligible; 4) the inclusion of rotation 
implies a downward shift of the whole distribution with an higher efficiency for the 
heavy s-only isotopes, leading to a flatter s-only distribution; 5) different prescriptions 
on convection or mass-loss produce nearly rigid shifts of the whole distribution. 

In summary, a variation of the standard paradigm of AGB nucleosynthesis would allow 
to reconcile models predictions with Solar System s-only abundances. Nonetheless, 
the EEPP cannot be definitely ruled out, because of the uncertainties still affecting 
stellar and Galactic chemical evolution models. 

Subject headings: Stars: AGB and post-AGB — Physical data and processes: Nuclear 
reactions, nucleosynthesis, abundances 
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1. Introduction 


Mass-losing Asymptotic Giant Branch (AGB) stars are the main souree of medium- 
and long-term gas returned to the interstellar medium (ISM). For this reason, they allow late 
episodes of stellar formation, thus prolonging the star-forming lifetime in many different Galactic 
environments. In addition, as a result of a eomplex eombination of internal nueleosynthesis and 
deep eonveetive mixing, the wind of AGB stars is heavily enriehed in both light (C, N, F, Na) 
and heavy elements. About half of the isotopes from Sr to Pb are produeed by AGB stars in t heir 


interior through a slow neutron eapture proeess ealled s-proeess (see, e.g.. 


BussoetaL 


1999 ). 


Moreover, the dust forming in their eool extended eireumstellar envelopes effieiently pollutes the 
ISM. Therefore, AGB stars play a fundamental role in the ehemieal evolution of galaxies. 

In this paper we diseuss the evolution of the heavy elements (A>90) in the solar 
neighborhood. Our main goal is to understand if the eurrent nueleosynthesis models provide a 
reliable evaluation of the ISM eontamination by AGB stars. The main and the strong eomponents 
of the s-proeess (A>90) are produeed by low-mass AGB stars, typieally 1.5<M/Mo < 3.0. Lighter 
s elements (A<90) are mainly synthesized by the s-proeess i n massive stars during eore He 


burning and shell C 


19921: 


turnin g (the so-ealled weak component; 


Pignatari et al. 


Kappeler et al. 


19891: 


Beer et al. 


20101) . Massive stars are also responsi 


eapture nueleosynthesis; for a review see 


Sneden et al. 


5le for the r proeess (rapid neutron 


20081). Most of the isotopes heavier than 


iron are produeed by both the s and the r proeess. However there exist a few isotopes that eannot 
reeeive any eontribution from the r proeess and, for this reason, are ealled s-only isotopes. An 
s-only isotope with atomie number Z is shielded by the r proeess due to the existenee of a stable 
isobar with Z-1 or Z-2. For this reason, the sequenee of /3 deeays that oeeurs at the end of the r 
proeess is interrupted before the s-only nueleus is reaehed. 


Galaetie Chemieal Evolution (hereinafter GCE) models obtained by eombining the s 
proeess eontribution of AGB stars (main and strong eomponent) and massive stars (weak s and 
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r process) have been studied by 


Travaglio et al. 


Travaglio et al. 


(11999 . 


2001 . 


2004h and 


Bisterzo et al. 


(l2014h . 


(l2004h firstly reported a deficit of the predicted Solar System abundances of Sr, Y 


and Zr (about -18%). These three elements belong to the first s-process peak in the Solar System 
composition, which corresponds to nuclei with magic neutron number N=50. After analyzing the 
possible uncertainties in their nucleosynthesis calculations, they concluded that this deficit would 
imply the existence of a missing s-process contribution, the so-called Light Element Primary 
Process (LEPP). Note that a different EEPP has also been invoked to explain the abundances of 


a large group of lig 


Montes et al. 


It elements with an important contribution from the r-process. Eor instance. 


(1200711 distinguished between ’’solar” and ’’stellar” EEPP, the latter found in 


metal-poor halo stars enriched by an r-process. Our findings are limited to the main s-process 
from AGB stars and, thus, we only focus on s-only isotopes in the solar nebula. Therefore, results 
presented in this paper do not provide any hint to certify (or exclude) the existence of a metal-poor 
primary EEPP, which could have equally well its roots in a sort of weak r-process^ _ 


he ne ed of an unknown pure s-process contribution has bee n also claimed 


(|2014t) basing on the analysis of the s-only isotopes (see also 


3y 


Kappeler et al 


Bisterzo et al 


201 ih . Indeed, 


in their chemical evolution models all the s-only isotopes with 90<A<130 are 


underestimated. As a matter of fact, the AGB yields used by 


Bisterzo et al. 


Travaglio et al. 


system atically 


( 20041) and 


1998 ) in which the 


20141) are based on post-process calculations (iGallino et al.l 
main neutron source (the pocket) is artificially introduced. The s process in low mass AGB 


stars is mainly due to the neutrons released by the ^^ C(a,n)^^Q reaction in a thin pocket 


o6r 


n. 


that forms after each third dredge up (TDU) episode dStraniero et al. 


1995 


Gallino et al. 


1998 


wStraniero et al. 


20061) . At present, a reliable evaluation of the extension in mass and of the 

ers 


(Herwie et al. 

1997; 

Denis senkov d 

Tou 

2003; 

Cristallo et al. 

the GCE models by 

Bisterzo et al. 

(2014 

), the extension of the 


task for AG 

B stellar model 

2009, 

2011; 

Eiu et al. 

2014) 


13 / 


fractions of and (the main neutron source and the main neutron poison, respectively) are 
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freely varied in order to reproduce the 100% of solar with an s-only distribution as flat as 
possible. These authors, however, did no t explo re the physical motivation at the base of those 


variations. More recently. 


Trippella et al. 


(12014) argued that, in stars with M<1.5 M©, magnetic 


fields are able to sha pe larger pockets than those characterizing more massive AGBs (see also 


Majorca et al 


201211 and suggest that this occurrence might have important consequences on the 


Solar System s-only distribution. However, their conclusions have to be verified with the support 
of a GCE model as well as that of evolutionary models that include the feedback of magnetic 
fields. 

The analysis of the LEPP problem presented in this work is based on a different approach. We 


verify if our EUll Network Stellar (PUNS, see 


Straniero et al 


20061 and references therein) yields. 


incorporated into a chemical evolution model for the solar neighborhood, can provide a reasonable 
fit to the Solar System s-only distribution. The adoption of a large nuclear network directly 
coupled to the physical evolution of the stars as well as our handling of the convective/radiative 
interface at the base of the convective envelope (i.e. where the forms) do not allow us to force 
our calculation to fit the absolute value of the Solar System s-only distribution. Notwithstanding, 
we can evaluated the effects on the AGB nucleosynthesis of different prescriptions for convective 
overshoot during the TDU, rotation-induced mixing, pre-AGB and AGB mass-loss rates and 
nuclear reactions efficiencies. 


The paper is structured as follow. We firstly describe our Galactic Chemical Evolution model 
and the stellar evolutionary code used to determine the proto-solar distribution for s-only isotopes. 
(Section [2] and Section [3l respectively). Then, in Section |4] we present our reference case, while 
in Sections [5] and we describe how AGB models uncertainties and GCE uncertainties affect our 
results, respectively. Our conclusions follow in Section |71 
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2. The Galactic Chemical Evolution Model 

We use a simplified GCE model for the solar neighborhood, defined as a eylinder of ~ 1 kpe 
radius at a distanee of ~8 kpe from the Galaetie eenter, adopting the standard formalism (e.g. 


Pagell120091) ■ Our GCE code is an update of that used to follow the evolution of light elements in 


previous studies (lAbia et al. 


19911 


19951) . The classical set of equations are solved numerically 


without the instantaneous recycling approximation (i.e. stellar lifetimes are taken into account) 
and assuming that at the star’s death its ejecta are thoroughly mixed instantaneously in the local 
ISM, which is then characterized by a unique composition at a given time. Thus, our predictions 
represent average values in time as this simplified approach cannot account for the scatter in any 
Galactic observable. Our main goal is to reproduce the absolute and relative isotopic abundances 
distribution of s-only nuclei at the Solar System formation (occurred 4.56 Gyr ago). These nuclei 
are ""Ge, ’^Se, '^’^^Sr, ^^Mo, '“"Ru, i°"Pd, “"Cd, “^Sn, i^s.noxe, 

i42Nd, I48.i50sjn, i 60 j)y^ ivoyb, i’Slu, '’^Hf, ^^^Os, 204p|, concentrate 

on those isotopes because they are produced only via the s-process and because an AGB origin 
is certain for those with atomic mass A > 96. We use the absolute isotopic abundances of the 
proto-solar nebula to normalize the output of our GCE model, which is stopped at that epoch. 
Note that those abundances differ from the current ones observed in the solar photosphere due to 
the impact of chemical settling. We obtain our proto-solar distribution by adopting the elemental 


abundances of 


Eodders et al 


procedure describe in 


(120091) and by computing a Standard Solar Model according to the 


Piersanti et al. 


(120071). 


The basic ingredients of the GCE model are described in the following. Eor the stellar yields 


see 


We adopt a standard Salpeter Initial Mass Eunction (IME), <1)(M) oc with A = 1.35 

in the mass range 0.1-100 Mq. Eor the Star Eormation Rate (SER) we have adopted the standard 
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Schmidt type law 'P(0 = acr^ (t), where CTgas is the surfaee gas density, a = 0.32 Gyr ^ and k = I 


We assume that the disk has been build up starting from an initial surfaee gas density (cTo) and by 
slow aeeretion of gas with primordial eomposition. Henee the initial abundanees of all the studied 
isotopes is set to zero. We adopt an exponentially deereasing gas aeeretion law f(t) oc ^“7^ 
infall preseription, eombined to the adopted SFR law, leads to a deereasing star formation history 
in the solar neighborhood. We set r = 7.5 Gyr sinee it has been shown that sueh a long timescale 


provides a satisfactory fi t to the observed ste 


solar neighborhood (e.g. 


Boissier & Prantzos 


ar Me tallicity Distribution (hereinafter MD) in the 


19991) . We have normalized t he infall rate f(t) by 


impo sing that the current observed total surface density is ~ 50 Mg pc ^ (see 


Goswami & Prantzos 


2000L and references therein for a detailed discussion). 


The main observables in the solar neighborhood which must be fitted are: 


• the current surface density of gas (13 + 3 Mgpc ; gas fraction 0.15-0.25%), stellar surface 
density (35 + 5 M© pc~^ ), total mass (~ 50 Mp, pc~^) as well as the current star formation 


rate (2-5 MgGyr ^pc ^) (iBoissier & Prantzos 


19991: 


Goswami & Prantzos 


20001 ) : 


the observed age-metallicity relation te.g. lCasagrande et al.ll201lh : 


the Type II and la supernova rates in the Galaxy (|Li 


[Fe/H] relationship in thin and thick disk stars (e.g. 


et al.ll201 ih. and the observed [O/Fel vs . 


Ramirez et al. 


20131: 


Nissen et al. 


2Q14t); 


• the observed MD o f long-lived G-type stars (ICasagrande et al.ll201llc lAdibekyan et al.ll2012 : 


Bensbv et al. 


2014 : 


the absolute a nd relative s-only iso topic abundances distribution at the formation of the 


Solar System (iLodders et al. 


2009). 


Recent studies (Roskar et al. 


2008 


Schonrich & Binnev 


2009 


Kubrvk et al. 


20131) have shown 


that the existence of gas and star migration across the disk of the Milky Way can significantly alter 
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the local observed age-metallieity relation and the stellar MD. One of the main eonelusions of 
those studies is that these observational eonstraints ean be properly interpreted only if migration of 
stars and gas is ineluded in GCE models. For instanee, reeent GCE models that inelude migration 
show that the average age-metallieity relation for stars locally born is generally flatter than the 
one ealeulated elassieally (i.e. without migration, as it is typieally done in one zone, ID GCE 
models). In particular, it implies that the Sun was not probably born locally, but it migrated from 
inner (more metal-rieh) Galaetie regions up to its eurrent position (r ~ 8 kpe). Stellar migration 
also introduees a dispersion in the observed abundanee ratios as a funetion of the metallieity (i.e. 
[X/Fe] vs. [Fe/H]). Although this dispersion seems to be generally small (cr < 0.15 dex), it might 
be larger for elements produced in low-mass, long-lived stars, like Fe or s-elements. Furthermore, 
the average gas metallieity in the ISM might differ from that of the loeal stellar population. The 
impaet of the gas and star migration in the observed s-elemen t distribution in the Solar System 


is beyond the seope of this study. We refer to speeifie studies (IKubryk et al. 


20131) for a detailed 


diseussion on the effeets of migration on the ehemieal evolution of the Galaxy. 

Note that we eannot a priori exclude that the simplifieations inherent our GCE eode mask some 
ehemieal features or introduee some biases in following the ehemieal evolution of s-only isotopes. 
Indeed, we are aware that more sophistieated models than the one zone GCE approximation 
adopted in the present work ean be eonstructed for the solar neighb orhood. We refer, for instan ce. 


to works ineluding the evolution of the h alo and of the thiek disk (IGoswami & Prantzos 


Kobavashi et al 


20001 : 


Mieali et al.l 


20001 : 


20131) . Different preseriptions for the SFR and infall/outflow 


of gas are typieally adopted for the evolution of these two Galaetie struetures, whieh are mainly 
eonstrained by their observed MD function. Nonetheless, we have ehecked that our results for 
the s-only isotopes abunda nee distribution are not affee ted when adding, for instance, the halo 


evolution (aeeording to the 


Goswami & Prantzos 


2000l prescriptions), provided that the initial 


metallieity for the disk evolution does not signifieantly exeeed [Fe/H]-1.0. Typieally, this is the 

maximum value of the metallieity for the halo reaehed in most GCE models after t~ 1 Gyr. 
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3. Stellar Models 


Stellar lifetimes, remnants masses and yields for low and intermediate mass stars (1.0< 
M/M© < 6.0) are derived from theoretie al evolutionary models eomputed with the FUNS 


evolutionary eode 


Straniero et al. 


2006 


. The stellar yields have been obtained by evolving 


models with different masses and initial ehemieal eomposition from the pre-Main Sequenee 
up to the AGB tip. In our models, the adopted AGB mass-loss rate has been ealibrated on the 
period-luminosity an d period-mass loss relations observed in Long Period Variable Stars (see 


Straniero et al. 


20061 and referenees therein). The atomie and moleeular opaeities in the eool 


envelope of AGBs aeeount for the variation of th e ehemieal composition as due to the occurrence 


of recurrent TDU episodes (iCristallo et al. 


200711 . During TDU episodes, the instability occurring 


at the inner border of the convective envelope is handled by adopting an exponential decay of 
the convective velocities. This makes the TDU deeper; moreover, as a by-product, we obtain 
the self-consistent formation of the po cket after each thermal pulse (TP) followed by TDU 


( Cristallo et al.l 


2009 


Straniero et al. 


2014) . The extension of the pocket varies from TP to TP 


following the shrinking of the region between the He and the H-shells (the so called He-intershell). 
The nuclear network used to follo w the physical and chemical evolution of our models has been 


presented in 


Cristallo et al. 


(1201 ih : it includes about 500 isotopes (from hydrogen to bismuth) 


linked by more than 1000 reactions. Such a network is directly included in the FUNS code, thus 
avoiding the use of post-process techniques. 

The main neutron source in AGB stars is represented by the ^^C(a,n)^^0 reaction, burning in 


'The FUNS code has been derived from the FRANEC code (Chieffi & Stranien 


Chieffi et al 


1998 ). 


19891: 
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- 1.5 -1 - 0.5 0 

[Fe/H] 


Fig. 1.— Selection of FRUITY yields for key s-process elements. Upper panel: (representative 

of the first s-process peak); intermediate panel: '^^La (representative of the second s-process peak); 
lower panel: (representative of the third s-process peak). See the on-line edition for a colored 

version of this figure. 
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radiative conditions at T~ 10^ K during the inter-pulse phases. An additional contribution conies 
from the activation of the ^^Ne(Q',n)^^Mg reaction at the base of the convective shells generated by 
TPs when the temperature exceeds 3 x 10^ K. While the former reaction dominates the s-proce^ 
nucleosynthesis in low mass AGB stars, the latter becomes fully efficient in stars with M> 3 
The models we use to calculate the AGB yields have di f ferent masses (1.0< M/M© < 6.0) and 


metallicities (-2.15<[Fe/H]<-i-0.15; 


Cristallo et al. 


2009L 


2011 


; Cristallo et al. in preparation). 


The corresponding yields are available on-line on our web repository FRUITTQ which represents 
our reference set. 

In Figure [U we report a selection of FRUITY net yield for some key s-process elements 


as representative of the first s-process peak, ^^^La as representative of the seconc 
andf^^Pb as representative of the third s-process peak). As already remarked in 


s-process peak 


Cristallo et al. 


(1201 ih . the largest yields are produced in the (1.5-3.0) M© mass range. Figure [U shows that low 
mass models (M<1.5 M©) marginally contribute to the global s-process production, since the 
TDU practically ceases to occur when the initial stellar mass drops below 1.2 M© (see also ^5.31) . 
Similarly, s-process yields from more massive AGBs (M>4.0 M©) are low, even if these stars 
may significantly contribute to the nucleosynthesis of some neutron rich isotope (as, for example, 
^’Rb and ^^Zr) due to the activation of the ^^Ne(Qr,n)^^Mg reaction. As expected, for stars with 
masses between 1.5 and 3.0 M©, the relative distribution of the three s-process peaks weakly 


^This mass limit depends on the metallicity. As a general rule, the minimum mass decreases 
with the metallicity. 

^http: //fruity, oa-teramo .inaf.it 
^A net yield is defined as 


X 


dM 

[Xik)-X\k)] - dt 

dt 


( 1 ) 


where dM/dt is the mass loss rate, while X(k) and X^{k) stand for the current and the initial mass 
fraction of the k-isotope, respectively. 
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depends on the mass, while it has a different behavior depending on the initial iron content. At 
large metallicities ([Fe/H]>-0.3), the s-process mainly populates the first peak (Sr-Y-Zr region). 
At intermediate metallicities, the second s-process peak (Ba-La-Ce-Nd region) presents its 
maximum. At low metallicities ([Fe/H]<-0.7) lead production dominates. 


In our GCE we adopt a simplified prescription by assuming that all stars with mass M 
> 8 M© explode as core collapse supernovae leaving behind a compact remnant such a neutron 
star of mass <1.4 M©, or a black hole in the case of most massive stars (M > 40 M©). In our 
calculations, we do not include the contribution of massive stars to the s-pro cess inventory. Those 


stars largely contribute to the production of s-only isotopes with A < 87 (see lPignatari et al. 


2010 . 


and references therein). Therefore, for those isotopes our prediction have to be considered as 
lower limits. Oxygen and iron yields from massive stars are instead needed in order to reproduce 


the average 


Nissen et al. 


Q/Fe1 vs. [Fe/H] relationship observed in unevolved sta rs (e.g. 


20141) . For that purpose, we use the yields published by 


Ramirez et al. 


Chiefh & Limongil (120041) . 


2013 


As far as it concerns the core-collapse supemovae contribution to the iron enrichment, we assume 
that on average each supernova ejects 0.1 of ^^Fe. On t he oth er hand, we adopt the type la 


supernovae explosion rate according to 


Greggio & Renz.inil (119831) in the framework of the Single 


Degenerate scenario for their progenitors. This corresponds to assume that a fraction of 2.5% of 
all binary system ever formed in the adequate mass range will provide an explosive outcome. This 
fraction value is set by fitting the observed current Galactic SN la rate and the [0/Fe] vs. [Fe/H] 
relationship. We also assume that, on average, f or each SNIa event an amount of ~ 0.7 M© of ^^Fe 


is ejected (e.g. 


Bravo & Martlnez-Pinedo 


2012 ). 


4. Reference case 

Our Reference case has been computed using the GCE model described in Oby adopting 
parameters values reported there. The model accounts for all the constraints mentioned above 
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Fig. 2.— Age metallicity relation (solid thiek eurve) eompared t o the average and +lo~ 


(dotted eurves) of t he observations in the solar neighborhood by 


Bensby et al. 


Casagrande et al 


(1201 it) and 


imits 


(l2014t) . The dashed curve refers to a GCE model computed with an increased SFR 


(+10% with respect to the solid thick curve). See Section!^ for details. 
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within the observational uneertainties. It is very well known that other reasonable ehoiees of the 
GCE model parameters (SFR law, IMF, ete.) might give similar results still in good agreement 
with the observational eonstraints. Due to its relevanee for our diseussion, we show in Fig. [2] 


the age-metallieity relation obtained in the Reference ease (thiek c ontinuous line). 


curves repres e nt the average and +lcr limits of the observations of 


Bensbv et al. 


Casagrande et al. 


he dotted 


1201 lh and 


(l2014h . Our model predicts a rapid increase of the ISM metallicity with time. 


reaching [Fe/H]« 0.0 at the epoch of Solar System formation and a continuous increase of [Fe/H] 
until now. 

In Fig. [2] we report the results of our Reference GCF calculation case, as obtained by 
using stellar yields included in the FRUITY database. Corresponding data are reported in Table 
[TJ In the upper panel we report absolute percentage s-only isotopic abundances obtained from 
our GCE model at the epoch of the Solar System formation. In this plot we did not add any 
contribution from the weak-s process. Thus, 100% means that an isotope is entirely synthesized 
by AGB st ars. For each i s otope , we also plot the corresponding solar abundance uncertainty 


reported in 


Fodders et al. 


(120091) . Dashed horizontal lines identify a +10% tolerance region 


representing th e current uncertaintie s in the estimated chemical abundances as due to nuclear 


cross sections (IKappeler et al. 


201 lh . In the lower panel of Fig. [3] we report the overproduction 


factors normalized with respect to ^^°Sm. In this case, unity means that an s-only isotope is 
over-produced (or under-produced) as ^^°Sm with respect to the corresponding solar abundance. 
The latter isotope has been chosen as refe rence since the entire s-process flux passes through it. 


1999j) . Also in this case, we highlight a 


making this isotope virtually unbranched (|Arlandini et al. 

±10% tolerance region. 

An inspection to Fig. [3] (upper panel) reveals an overall overproduction of s-only isotopes with 
A>96 (~ 145%), more evident in the region 128<A<204. Thus, on a relative scale, lighter s-only 
isotopes are underproduced with respect to the heaviest ones (see lower panel of Figure [H). 
However, our relative distribution can be considered flat if current uncertainties (observational and 
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50 100 150 200 


A (Atomic mass) 


Fig. 3.— Reference GCE calculation case. Upper panel: percentage s-only isotop k abundances 


at the epoch of the Solar System formation. Solar abundance errors are taken from 


Fodders et al. 


(120091) and reported in Table [U (error bars smaller than symbols are not reported). Lower panel: 
overproduction factors normalized to the ^^°Sm one. Dashed horizontal lines identify a +10% 
tolerance region in both panels. See text for details. 
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nuclear) are taken into account. Si nce we do not assume any ad hoc re-scaling of the pocket 


at odds with 


Travaglio et al. 


(1200411 and 


Bisterzo et al. 


(1201411 who claimed a missing contribution 


to light s-only isotopes, we obtain super-solar percentages for all s-only isotopes with a sure AGB 
origin (A>96). Thus, in the following we investigate if there is the possibility to decrease the 
overall Galactic s-only production and if a larger depletion efficiency can be found for the heavier 
s-only isotopes (128<A<204). This exploration is carried out in the next Section by studying 
current uncertainties affecting stellar models. 

As starting point of our analysis, however, we want to verify if GCE models confirm that the 
bulk of the s-process comes from AGB stars with masses 1.5 <M/Mq <3.0 (see previous Section). 
The contribution from AGBs with M<1.5 M© will be analyzed in ^5.31 In order to quantify the 
contribution to the Solar System s-only distribution from Inter mediate Mass Stars AGBs (i. e. stars 


with initial masses M> 4 M©, hereinafter IMS-AGBs; see also 


Karakas & Lattanzio 


2014h . we 


run a GCE model by setting to zero the yields of those objects (hereinafter No IMS case). Results 
are shown in Eigure|4]and reported in Table [B On average, the IMS contribution to the Solar 
System s-only distribution i s marginal (on averag e 6%). Thus, even if our IMS-AGBs present 


tiny '^C-pocket after TDUs dStraniero et al. 


1201411 . their contribution, once weighted on the IME, 
is small. Eor the lightest s-only isotopes (form ^°Ge to *^Sr), the relative IMSs contribution is 
larger, due to the more efficient activation of the ^^Ne(a,n)^ ^Mg source. Note, however, that those 


isotop es are mainly synthesized by the weak-s co mponent (IKaeppeler 


201011 . Thus, we basically confirm the finding of 


et al. 


1994 


Pignatari et al. 


Bisterzo et al.! (20141 that intermediate mass 


AGBs marginally contribute to the Galactic chemical evolution of s-only isotopes. 


5. Stellar Models Uncertainties 

Despite great strides made by stellar modelers in the last two decades, our understanding of 
the AGB phase is still hampered by large uncertainties. A few physical details can be constrained 
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Fig. 4.— As in Figure[3l but including a GCE calculation without the contribution from IMS AGBs 
(crosses). The Reference ease (open dots) is shown by eomparison. We have omitted error bars for 
elarity. 
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by theory and, therefore, the adoption of phenomenological models is often the only way to 
describe a specific physical process. Thus, it is not surprising that a large part of the extant 
models still include a set of parameters sometimes rather free, sometimes (partially) constrained 
by observations. Despite all these limitations, we try to evaluate the effects that current stellar 
modelling uncertainties have on s-process yields. In the present work we focus on some key 
physical processes (such as rotation, convection and mass-loss) and on the efficiency of some 
nuclear processes. 

Due to the large number of models included in the FRUITY database, we compute a reduced 
number of (M,Z) combinations by analyzing once a time each of the above mentioned physical 
processes and we derive corrective factors to be applied to models with similar masses (M) and 
metallicities (Z). Such a procedure does not introduce biases in our conclusions because we focus 
our attention on those (M,Z) combinations where, according to our previous experience, major 
effects are expected. 


5.1. Rotation 


FRUITY AGB stellar models are representative of the intrinsic carbon stars observed in the 
disk and in the halo of the Milky Way. However, a comparison between our theoretical curves and 
spectroscopic data show s that, at fixed metalli city, our models do not cover the observed spread 


in the s-process indexes. 


Piersanti et al. 


(1201311 recently demonstrated that a variation in the initial 


Zero Age Main Sequence (ZAMS) rotational velocity {'^zams^ determines a consistent spread in 
the final surface s-process enhancements and spectroscopic indexes in stars with the same initial 
mass and metallicity. Rotation-induced instabilities (in particular the Goldreich-Schubert-Fricke 
instability and meridional circulations) modify the mass extension of both the and the 
pockets and their relative overlap. This is shown in Figure[5l where we report the and the 
mass fractions in the upper layers of the He-intershell after the 4'^ TDU episode of a 2 Mq model 
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with Z=10 ^ ([Fe/H]=-0.15) and km/s. We plot chemical profiles at the end of the 


formation of the pocket (dotted lines) and at the beginning of the neutrons release by the 
'^C(Qf,n)^^0 reaction (solid lines). The abundance of is also plotted to testify the starting of 
neutron capture processes. With respect to a non-rotating model the average neut ron-to-seed rati o 


Piersanti et al. 


decre ases and the production of s-process elements is lower (see also Figure 3 in 
2013(1 . This is a consequence of the higher abundance of a very strong neutron poison, in 


the pocket. This also implies that light-s elements are less depleted than the heavier ones. It 
is worth mentioning that the inclusion of rotation does not s ubstantially affect the efficiency of 


Piersanti et al 


2013 or the FRUITY 


TDU, as testified by the almost unaltered surface [C/Fe] (see 
database). 

Rotational velocities of Main Sequence stars of spectral classes A and F span on a quite large 
range and they can be as high as 300 km/s. On the other hand, asteroseismol ogy measurements 
seem to indicate that the cores of Red Giant stars rotate quite slowly (see e.g 


Mosser et al. 


201 


This discrepancy is normally attributed to a particularly efficient transfer of angular momentum 
from the inner zones to the convective envelope or to magnetic braking. In our models we do not 
account for such an effect, but to compensate it we use low ZAMS rotation velocities. Thus, we 


assume km/s for models with M<2.0 M© and a slightly larger value for models with 

'ZAMS km/s). Due to the marginal contribution to the bulk of the s-process 


2.0<M/Mo <4.0 (vf -- 


from IMSs (see Ml), we do not apply rotating corrective factors to the yields of more massive 
AGBs (M> 4 Mq). In Figure [6] we compare our Reference case with a GCE calculation based 
on stellar models including the effects of rotation (hereinafter Rototzon case). The corresponding 
data are reported in Table [IJ As expected, we find a general decrease of the absolute s-process 
abundances, the depletion factors increasing for larger atomic masses. With the exception of 
^^^Ba, s-only nuclei show absolute sub-solar percentages. The lightest s-only isotopes (up to ^®Sr) 
are less depleted than the heavier ones. In fact the reduced neutron exposure (due to the partial 
overlap between the and the pockets) leads to the synthesis of isotopes closer to the iron 
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U/U^ 


Fig. 5.— and profiles after the 4'^ TDU in a 2 Mq rotating model with Z=10^^ ([Fe/H]=- 

0.15) and y^zAMs~^^ km/s (solid eurves). The profile (multiplied by a faetor 1000) is also 
displayed. We show and profiles at the end of the formation of the poeket (dotted 
curves) and when neutrons start being released (solid curves). See the on-line edition for a colored 
version of this figure. 
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Fig. 6.— As in Figure IH but including a GCE model with rotating AGB models (crosses). 
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seeds (^^Fe), to the detriment of the heavier ones. This is even more evident when looking to the 
relative overproduction factors (lower panel of Figure!^. On a relative scale, light s-only isotopes 
gain more than a factor 2 with respect to the Reference case, while those in the atomic mass range 
96< A <124 are now within (or even above) the tolerance region. 

Obviously, a different choice of the initial rotational velocities would lead to a different s-only 
distribution in both the absolute and relative scales. Thus, in principle a better fit could be found. 
However, due to the other uncertainty sources affecting ste _ 
to the treatment of rotation in ID evolutionary codes (see 


lar models, in particular those related 

Piers anti et al. 

2013 

and references 


therein), we prefer to highlight general effects related to a physical input (as rotation) more than 
to provide a detailed specific recipe to obtain the desired fit. 


1968 


5.2. Convection 


n our models, according to the prescriptions of the Mixing Length Theory (MLT; 


Cox 


, convective velocities are proportional to the difference between the radiative and the 


adiabatic temperature gradients. Thus, in presence of a smooth adiabatic temperature gradient 
profile, the convective velocity is 0 at radiative/convective interfaces. However, when the 
H-rich envelope penetrates in the He-rich region, which is characterized by a lower opacity, a 
non-zero convective velocity is found at the inner border of the convective envelope^. This is the 
standard picture of a TDU episode. As a consequen ce of this abrupt chan ge in the opacity, the 


radiative/convective interface becomes unstable (see 


Straniero et al. 


20061 for details). However, 


the steep pressure gradient should limit the penetration of such an instability and, thus, the average 
convective velocity should rapidly drop to 0. We mimic this behavior by assuming that convective 
velocities follow an exponential decay law below the convective envelope. This has two major 


^We remind that the radiative gradient is proportional to the opacity. 
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conse quences: the TDU episode is deeper and, later, a '^C pocket develops (see iCristallo et al. 


20091 for a detailed discussion and for a comparison with techniques used by other groups to 
handle the formation of the pocket). 

In our FRUITY models, the penetration of protons is inhibited below 2 Hp from the formal 
Schwarzschild Boundary (hereinafter SB). In order to explore the sensitivity of stellar yield 
and, hence, of GCE calculations on such an assumption on mixing efficiency, we computed the 
same FRUITY non-rotating 2 M© Z=10^^ ([Fe/H]=-0.15) model, but allowing the partial mixing 
below the SB down to the layer where the convective velocity is 10“^^ times the value at the SBg 
(hereinafter Tail case). In the upper panel of Figure 0 we plot the and abundances after 
the 3'“^ TDU episode of the FRUITY model (solid curves) and the Tail model (dotted curves). As 
it can be easily derived, the integrated mixed below the SB is almost the same for the two 
cases (7.6xl0~^ M© vs 7.8xl0~^ M© for the FRUITY and Tail models, respectively), while the 
integrated is 50% larger in the Tail model with respect to the FRUITY one (2.3x10“^ M© vs 


3.4x10 ^ M©). This means tha t the effective (i.e. the that effectively contributes to the 


s-process; 


13 




Cristallo et al. 


201 ih is nearly twice in the Tail model (from 9.2x10 ® M© to 1.8x10 


M©). As a consequence, the overall s-process production increases, as testified by the curve in the 
lower panel of Figure |71 where we plot the differences in the final surface enhancements between 
the Tail and the FRUITY models. Corresponding data are reported in Table [B As expected, 
light elements (Z< 28) are not affected by the changes in the tail profile, while the three 
s-process peaks show larger surface enrichments (about -1-30% for Is and hs and -1-60% for lead). 
This additional contribution comes from the portion of the extended pocket characterized 
by 5 X 10 ^ < X(^^C) < 1 X 10 ^ (see Figure [T]). In our previous models, such a contribution is 
suppressed since it lies in correspondence to the drop of the profile. Interestingly, elements 
normally associated to the weak component (Ge-Ga) result strongly enhanced with respect to the 


®This roughly corresponds to (2.2-2.4) Hp. 
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M/Mg 



Fig. 7.— Upper panel: and profiles after the 3''“^ TDU of a 2 Mq, Z=10“^ ([Fe/H]=-0.15) 

models with different preseriptions for the radiative/eonveetive interfaee treatment. Lower panel: 
elemental surfaee differences between the two cases shown in the upper panel. See the on-line 
edition for a colored version of this figure. 
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O 



M/M^ 


Fig. 8.— Selected key isotope profiles in the pocket layers after the 3''“^ TDU of a 2 Mq and 
Z=10“^ ([Fe/H]=-0.15) model. Thick and thin curves refer to the Tail and FRUITY cases, respec¬ 
tively. See the on-line edition for a colored version of this figure. 
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Fig. 9.— As in Figure |4l but including a GCE model based on AGB models handling in a different 
way the radiative/convective interface at the base of the convective envelope (crosses). See text for 
details. 
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FRUITY model (up to 60%). This is due to the contribution from the inner tail of the pocket, 
where neutron densities are lower and, thus, less massive isotopes are synthesized (see, e.g., the 
^^Ga profile in Figure [H]). In the Tail model, neutron rich isotopes (as ^®Zr), normally bypassed by 
the s-process main path, are not enhanced, but even mildly depleted with respect to the FRUITY 
model. This is a consequence of the larger mass extension of the pocket of the Tail case. 
During the radiative burning, in fact, neutron-rich isotopes are destroyed more than produced 
(see, e.g., the ^^Zr profile in FigurelH). 

In the upper panel of Figure |9] we show GCE absolute percentage s-only abundances obtained 
by using AGB models with the new prescription for the lower boundary of the pocket {Tail 
case). For s-only isotopes with A > 96 we find an increase of s-process absolute percentages (on 
average -1-30%), with similar enhancements for light and heavy s-only isotopes (as testified by the 
similar relative distributions reported in the lower panel). Lighter s-only isotopes (A< 87) result 
more enhanced with respect to the heavier ones due to the contribution from the inner tail of the 
pocket. In summary, we find that larger -pockets do not strongly modify the shape of th e 


s-only distribution, but sizeably affect their absolute values (see also Figure 4 in 


Bisterzo et al. 


2014t) . We also find that the Tail case is able to nearly 


reproduce the entire Galactic production 


of ^^Sr and '^'Sr, in agreement with previous findings (iTrippella et al. 


2014h . However, when 


compared to other s-only nuclei, Sr s-only isotopes are underproduced and, therefore, in our GCE 
model a certain contribution from the weak-s process is still needed. 

Note that a rigid shift (in both directions) could also be obtained by assuming a different a 


parameter of the MET. Sue 
Standard Solar Model (see 


1 a parameter is cali brated by reproducing the solar properties with a 


Piersanti et al. 


2007 


for details). However, there is n o specific reason 


to adopt the same a for all t he stellar evolutionary phases (see the discussion in 


2014i . 


Straniero et al. 


Cristallo et al.l (120091) already showed that a reduction of the MET parameter in AGB stars 


leads to a decrease of the s-process yields. Since this variation does not depend on the metallicity, 
we expect a corresponding rigid shift in the output of a GCE. 
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5.3. Mass-loss 

The poor theoretical knowledge of the stellar mass loss history represents one of the main 
uncertainties in the computation of AGB stellar models. Low and intermediate mass stars lose the 
majority of their mass during the Red Giant Branch (RGB) and the AGB phases. 


In ID stellar evolutionary codes, the mass-loss rate du ring the 


parameterized according to the formulation proposed by 


Reimers 


^GB p hase is commonly 

J1975I1 : 


M = 4x 10^'" — 
gR 


( 2 ) 


where M is in units of Mo/yr and other quantities are in solar units. The uncertainty affecting 
this formula was orig i nally quoted by Reimers to be at least a factor 2 either way. Later, 


Fusi-Pecci & Renzinil (Il976h introduced a normalization constant in order to reproduce the 
Horizontal Branch (HB) morphology of Globular Clusters (/7 k=0.4j)- Depending on the mass lost 
during the RGB phase (and thus on the value of ? 7 r), stars attain the AGB phase with different 
envelope masses. Thus, in principle, the RGB mass-loss could have an effect on the subsequent 
AGB nucleosynthesis. Those effects are expected to be important for low mass stars (M<L5 
Mq) because they spend more time on the RGB phase with respect to larger masses. Moreover, 
their envelopes are thinner with respect to more massive stars and, therefore, even a small amount 
of material lost (e.g. 0.1 M©) c an produce sizeable e ffects on the occurrence of TDU in the 


subsequent AGB phase (see e.g. 


Straniero et al. 


20031). To properly determine the effects of RGB 


mass-loss rate on AGB nucleosynthesis (and thus on the solar s-only distribution), we calculate 
a set of M=L3 Mq models at different metallicities with rjR = 0.2. In Figure ITOl we report the 
variations of the surface abundances (A,) with respect to the corresponding FRUITY cases. We 
find that A, are larger at low metallicities (in particular for the heaviest s-only isotopes). This is 
due to the fact that at large metallicities this mass experiences a few TDU episodes even using a 


^FRUITY models adopt this value. 


















-29- 


milder RGB mass loss rate. Thus, the final surface s-process enhancement is, in any case, low. By 
comparison, we also report data relative to a M= 2 Mq, Z=10~^ ([Fe/H]=-0.15) model. The low 
variations found in this case confirm that for massive enough AGB stars a reduced mass-loss rate 
during the RGB phase has practically no effect. In Figured!] we report a GCE model computed 
with rjn = 0.2 in stars with M< l.SM© (hereinafter Reimers case). We find minor variations in the 
s-only distribution (see also Tabled]), with slightly larger enhancements for the heaviest s-only 
isotopes (A > 128). Our results reinforce the evidence that the major contributors to the Solar 
System s-process inventory are AGB stars in the mass range (1.5-3.0) M©, as already inferred in 
Section [3] Their nucleosynthesis is strongly affected by the rate at which they lose mass during the 
AGB. A viable method to estimate AGB mass loss is based on the observed correlation with the 


pulsation period (iVassiliadis & Wooi 


19931) . Since the evolution of the pulsation period depends 


on the variations of radius, luminosity and mass, this relation provides a simple method to estimate 
the evolution of the mass loss rate from basic stellar parameters. In our m odels, the AGB mass 


OSS is determined according to a procedure similar to the one adopted by 


Vassiliadis & Wood 


19931), but revising the mass loss-period and the period-luminosity re 


more recent infrared observations of solar metallicity AGB stars (see 


ations, taking in to account 


Straniero et al. 


2006 and 


references ther ein). It has been demonstrated that AGB mass 


the metallicity (IGroenewegen et al. 


2007 


Lagadec et al. 


OSS rates are mildly dependent on 


2008h and, thus, we applied the same 


period-mass loss relation for all AGB models present in the FRUITY database. Notwithstanding, 
it is worth to note that, when a fixed period is defined, observational data show a quite large 
scatter. In a period-mass loss plot, a theoretical curve constructed reducing by a facto r 2 the 


mass -loss rate at a fixed period still lays within the observed spread (see Figure 8.10 of 


Cristallo 


2006h . This still holds for a mass-loss rate increased by a factor 2. In order to quantify the effects 


on the s-only distribution induced by a variation of the AGB mass-loss rate, we compute some 
AGB models with a milder and stronger period-mass loss relations. In Figure [l2] we show the 
results on the final surface distributions of 2 M© stellar models at various metallicities. In the 
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Fig. 10.— Differences in the surface chemical distributions of 1.3 Mq stars at various metallicities 
computed with a reduced RGB mass loss rate (REI: ? 7 « = 0.2) with respect to the corresponding 
FRUITY models (FRUITY: rjR = 0.4). A 2.0 M© star is also reported by comparison. See the 
on-line edition for a colored version of this figure. 
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Fig. 11.— As in Figure IH but including a GCE model with a reduced RGB mass-loss rate in low 
mass stars (crosses). 
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Fig. 12.— Differences with respect to FRUITY models in the final surface chemical distributions 


of 2.0 Mo stars at various metallicities computed with an i ncreased or decreasec 


rate. Our standard AGB mass-loss rate has been described in 


Straniero et al. 


2006 


AGB mass loss 
Negative differ¬ 


ences are obtained with an increased mass-loss rate (specularly, positive differences are found for 
models with a reduced AGB mass-loss rate). As expected, elements whose AGB production is neg¬ 
ligible show null differences (i.e. they have the same final surface abundances of the corresponding 
FRUITY models). See the on-line edition for a colored version of this figure. 
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Fig. 13.— As in Figure 01 but including a GCE model obtained by assuming a lower AGB mass 
loss rate (crosses). See text for details. 
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plot. A, represents the difference between models computed with the standard M-period relation 


(Straniero et al. 


2006h and the modified ones. Obviously, positive differences are obtained with a 


milder mass loss rate, while negative differences with the stronger one. Heavy elements surface 
variations are below 0.1 dex (25 %) for the whole s-process distribution, being slightly larger at 
low metallicity. Thus, we expect that a modified M-Period relation in the AGB phase produce an 
almost rigid shift (upward or downward, depending on the adopted mass loss law) of the s-process 
isotopes. In order to verify this statement, we compute a GCE model with a milder M-Period 
relation during the AGB phase (hereinafter Mloss AGB case). Results are shown in Figure [T3l 
corresponding data are reported in Table [B As expected, for s-only isotopes with A>96 there is 
an almost rigid upper shift of solar percentages (~25%). In summary, a rigid shift (upward or 
downward) of the s-process isotopic inventory can be obtained by adopting a different prescription 
for the AGB mass loss rate within the intrinsic observed scatter in the M-period relation. 


5.4. Efficiency of nuclear processes 


In previous Sections we demonstrated that different prescriptions on physical processes can 
lead to appreciable variations of the s-only inventory. In th is Section we concentrate on strong and 


weak nuclear processes. We refer to 


Cristallo et al. 


(1201 ih for a list of the adopted reaction rates 


in FRUITY models. Here, we focus on the uncertainties affecting the rates of: 


• nuclear processes determining the abundances of s-only isotopes close to s-process 
branchings; 

• neutron sources in AGB stars, i.e. the '^C(Q',n)^^0 and the ^^Ne(a,n)^^Mg reactions; 

• the major neutron poison in AGB stars, i.e. the ^"^N(n,p)^^C reaction. 








-35- 


z 



Eu 

Sm 




► 

A 


Fig. 14.— s-process main path in the region of and branching points. 

By means of the first test we can quantify local variations of s-only isotopes, while the others 
allow us to determine if nuclear processes are able to shape the whole s-only distribution. 


5.4.1. s-process branchings 

We focus on the branchings at ^^"^Cs and at ^^^Eu, which determine the surface abundances 
of and ^^^Gd (overproduced and underproduced with respect to ^^°Sm in our GCE models, 
respectively). In Eigure [T4l we report the main s-process path in the regions of the nuclide 
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chart corresponding to the two s-proeess branehing points. The unstable isobars have p deeay 
timeseales of the order of years (2.1 yr and 8.8 yr in laboratory eonditions for ^^^Cs and 
respeetively). Thus, their deeays are faster than eorresponding neutron eapture during the radiative 
burning, but long enough to allow the opening of s-proeess branchings during the convective 
^^ Ne(cK,n)^^Mg burn ing. Direet measurements of the ^^^Cs(n, 7 )^^^Cs reaetion is prohibitive 


(Patronis et al 


20041). while for the neutron eapture on ^^^Eu only a dated activation measurement 


by adopting the uneertainties reeently provided by 

Rausehei 

(2012 

) (+10% anc 

±50% for i34Ba 

and ^^"^Gd, respeetively). The (5 deeays rates are taken 

fom 1 

’akahashi & Yokoi 

(1987 

), while the 


corresponding uneertainties (a faetor 3) are taken from 


Gorielv ( 


1999 ). 


In Eigure[l5]we report the differences (filled dark eircles) in the yields of a 2 M© Z=10“^ 
([Ee/H]=-0.15) model with respect to a ERUITY model by modifying neutron cross sections and 
yS decays in the following way: 


• upper limit of the ^^"^Cs neutron eapture eross seetion; 

• lower limit for the deeay rate of ^^“^Cs into 

• lower limit of the ^^^Eu neutron eapture cross section; 

• upper limit for the decay rate of ^^^Eu into ^^^Gd . 


These choices aim at minimizing the ^^^Ba produetion and at maximizing the ^^^Gd production. 
The plotted quantities are normalized to variations in ^^°Sm yields (thus, unity means no variation 
with respeet to ^^°Sm). We find a 5% reduetion of '^"^Ba yield and a 30% inerease of ^^"^Gd yield. 
Note that these numbers refer to 2 M© models: for more massive AGB stars (e.g. 3.0-4.0 M©) 

^This corresponds to adopt the upper limit of the ^^^Cs lifetime, whieh is defined as the inverse 


of d (and thus of the rate). 



































1.4 


M = 2 Mq Z=1x 10-2 ([Fe/H] = -0.15) 



Fig. 15.— Differences in the yields of a 2 Mq model with Z=10“^ ([Fe/H]=-0.15) obtained by 
varying strong and weak reaction rates (dots: variations of nuclear processes efficiencies in cor¬ 
respondence of ^^"^Cs and branching; triangles: variations of neutron source cross sections; 
squares: variations of neutron poisons cross sections). Differences are normalized to variations in 

150 


Sm yields. See text for details. 





Overproduction factor to '“Sm % contribution by AGBs at t 


-38- 


,0 103 F 


10° 


T- 1 -1- 1- 1-1- 1 -1- 1 - 1 -1- 1 -1-T 


9s» 


0 


103 = = = = = = = z f zzzzzz f z®z z z-4 


• » 


9Q 


10' r 

■ Nuclear 


I I I 


I I I 


50 


100 


150 


200 



Fig. 16.— As in Figure IH but including a GCE model with different prescriptions on selected 
strong and weak nuclear processes (crosses). See text for details. 
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these effects are larger. 


5.4.2. s-process neutron sources 

Then, we verify if the solar s-only distribution is modified when adopting recently published 
rates for the two major neutron sources in AGB stars, i.e. the ^^C(Qr,n)^®0 and the ^^Ne(Qr,n)^^Mg 
reactions. While the first reaction releases neutrons in radiative conditions during interpulse 


13 / 


^l6r 


we used the value proposed by 

La 

Cognataetal. (2013). while for the ^^Nelor.nl^^Mg the 

value suggested bv 

Loneland et al. 

(2012 

) is adopted. With respect to our reference rates 

(Drotleff et al. 1993 

and 

Kaeppeler et al. 

1994. respectivelvl. both of them are about 20% 


higher at the temperatures of interest. The combined effect induced by the new '^C(a,n)^®0 
and ^^Ne(a;,n)^^Mg reactions is an overall slight increase of the whole s-only distribution. 
This derives from the fact that with an higher '^C(a,n)^®0 reaction rate the fully burns in 
radiative conditions, while when using the reference rate s ome of the in the first pockets 


can be engulfed in the convective shells generated by TPs (ICristallo et al. 


2009D- When is 


engulfed and burns convectively, only neutron-rich isotopes as ^‘^Fe and ^^Zr are synthesized. 
This is confirmed by the strongly reduced ^^Zr abundance (open pentagon) we obtain in this 
models despite the increased ^^Ne(Q:,n)^^Mg reaction rate. On a relative scale, we notice a 
marginal decrease of the lighter s-only isotopes with respect to the heavy ones (triangles in Figure 


[T5l) . The increase o f the ^^Ne 


production (see also 


Liu et al. 


Q-^n) 


25 


Mg reaction does not produce sizeable effects on the '^"^Ba 


2014) . which is at the same level of s-only isotopes with 96<A<130. 


Such a result further confirms that this reaction is only marginally activated in low mass AGB stars. 


^A similar decrease is also found for ®°Fe. 
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5.4.3. s-process poisons 


Major neutron poisons in AGB stars are the and the ^^Al(n,p)^^Mg reaetions, 

working in poekets and during TPs, respeetively. In eonsideration of the weak aetivation of 
the ^^Ne(a,n) ^^Mg neutron source, we c oneentrate on the first reaction only. Our reference rate 


is taken from 


Koehler & O’brienI (Il989h . In Figure [T5] we report the variations in s-only isotopes 


yields (squares) by considering an increased rate of 10%. We find a general decrease of heavy 
s-only isotopes, which translates in a general overproduction of light s-only isotopes with respect 
to (-1-8% on average). This is due to the fact that an increased poison effect reduces the 
s-process efficiency and, thus, its capability to by-pass the bottleneck at Z=50. Specular results 
are expected when considering the lower limit of the '^N(n,p)^^C reaction. Note that this effect is 
less relevant for higher masses (in which the ^^Ne(Q',n)^^Mg is more efficiently activated) and it 
practically vanishes at low metallicities (where the Z=50 bottleneck is more easily by-passed due 
to the larger neutrons-to-seeds ratio). 


5.4.4. Effects on a GCE model 

In Figure [T^ we report the results of a GCE model in which we take into account for all of 
the afore-described modified rates (hereinafter Aiic/ear case). Corresponding data are reported in 
Tabled! Variations in the cross sections of neutron sources and of the major neutron poison in 
AGB stars do not lead to significative changes in the global s-only isotopes distribution. However, 
this is not the case for the solar abundances of s-only isotopes close to s-process branchings. 

In fact, we find that the solar ^^^Ba and ^^^Gd percentages decrease and increase by more than 
20%, respectively. However, we remark that the uncertainties in the [3 decay rates (which mainly 





-41 - 


determine the differenees showed in Figures [l5] and [H]) are rough estimates and, thus, larger 
isotopie variations eannot be a priori exeluded. A further theoretieal nuelear analysis on this topie 
would be highly desirable. 


6. Galactic Chemical Evolution Models Uncertainties 


It is important to remind that also GCE models strongly depend on the adopted inputs, sueh 
as the SFR, the IMF or the type la supernovae evolutionary seenario (see SeetionlH). In faet, eaeh 
of these quantities influences the amount of metals locked or released by stars at different epochs. 
The relevance of their impact would depend on how much their variations (within the current 
uncertainties) would affect the derived age-metallicity relation. For instance, a faster increase of 
the ISM metallicity would imply a lower contribution from metal-poor stars, because there would 
be less time to form them. Thus, the contribution from metal-rich AGB stars to the Solar System 
s-process distribution would be larger and, consequently, t he produ ction of s-only isotop es with 


96 < A < 124 would result increased (Majorca et al. 


2012 


; see also 


Trippella et al. 


2014). In fact. 


the higher the iron seeds number the lower the atomic mass of the synthesized s-process nuclei. 
We plan to systematically study the impact of different choices of the GCE input parameters on 
the s-only distribution in a forthcoming paper. Hereafter we only show the effects that a variation 
of the SER has on the solar s-only isotopic distribution. 

Observations of various SER indicators in galaxies reveal that star formation occurs in different 
ways, depending on the galaxy type. There is no theory to predict star formation on large scales 
in a galaxy given the many physical ingredients that may affect the SER. In Eigure [2] we show 
the age-metallicity relation obtained by assuming an increased SER at all epochs (-1-10%; dashed 
curve) with respect to our Reference case. Here an increased SER mimics a higher [Ee/H]. 
Corresponding data are reported in Table [B We will refer to this case as SER-l-10. Actually, we 
have just changed the a parameter in the Schmidt’s law by 10%, because with this choice it is still 
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Fig. 17.— As in FigurelH but including a GCE model with an increased SFR of +10% (crosses). 
See text for details. 
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possible to account, within observational uncertainties, for all the solar neighborhood observables 
mentioned in SeetionlH A larger variation of the SFR would imply a new calibration of the GCE 
model itself. In this ease, however, it would be diffieult to disentangle the effeets related to the 
ehange in the SFR from those conneeted to the new parameter set adopted to fit again observables. 
As shown in Figure [JTl the variation of the SFR has an appreeiable effect on the s-only isotopes 
distribution. We notiee a slight increase of light s-only isotopes and a more consistent deerease of 
the heavy ones. As a eonsequence, on a relative scale light s-only isotopes with 96 < A < 136 are 
overproduced with respeet to the heavier ones by 18% on average. 


7. Discussion and Conclusions 


In this paper we verify if our FUNS stellar yields (available on the FRUITY database), 
used in a Galactic chemical evolution model, ean reproduce the distribution of s-only isotopes 
eharaeterizing the proto-solar nebula. Those nuelei are only synthesized by the s-proeess and, 
thus, are exeeptional markers of the evolution of pa st Galaetie AGB populations. At odd with 


previous studies based on post-process calculations (IBisterzo et al. 


2014 


Trippella et al. 


2 ^, 


we use in our analysis AGB stellar yields obtained by means of stellar evolutionary calculations 
fully coupled to an extended nuelear network. 

In our GCE models, we find that the contribution to the Solar System s-only distribution from 
low mass AGB stars (M< 1.5 Mq) as well as from intermediate mass AGB stars (M> 4 Mq) is 
marginal. Thus, we confirm that the bulk of the s-process eomes from AGB stars with masses 
(1.5 - 3.0) Mq. Another major result of this study is that, within the eombined uneertainties, we 


do not miss any contribution t o the Solar Sys t em s-o nly c 


96 < A < 124, as elaimed by 


Travaglio et al. 


(12004 and 


istribution in f 


Bisterzo et al. 


le ato mic mass range 


2014 ). Our referenee 


GCE model, in fact, predicts an overall super-solar s-only distribution (—1-45% on average). 
When observational and nuclear errors are taken into aeeount, the distribution relative to '^°Sm 
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can be considered flat, even if a lower production is found for s-only nuclei with 96 < A < 124. 
We investigate if current uncertainties affecting stellar models can lead to a better fit to the Solar 
System s-only isotopic distribution. The inclusion of rotation in our stellar models implies a 
general suppression of the s-process, with larger depletion factors for the heaviest s-only isotopes. 
On a relative scale, this implies a larger contribution to light s-only isotopes and, thus, a flatter 
s-only distribution. Different prescriptions for convection efficiency and for the treatment of the 
unstable inner border of the convective envelope during TDU episodes produce nearly rigid shifts 
of the entire s-only distribution. The same result can be achieved by adopting a different mass-loss 
rate during the AGB phase. Current nuclear uncertainties affecting strong and weak reactions 
allow for important improvements in the determination of some s-only isoto pes (as ^^^Ba and 


154 


Gd) . The need of revised /3 decay rates with respect to those published by 


Takahashi & Yokoi 


1987h is highly compelling. 


In the past, the nucleosynthesis of s-only isotopes has been closely related to that of ^°^Pb. 
Although such a nucleus is not a pure s-process isotope, a large per centage of its solar abundance 


is ascribed to the s-p rocess, the estimates varying from about 85% (ICowan et al 


(|Bisterzo et al. 


1999) to 98% 


20141) . Our reference model slightly overestimates its absolute solar abundance 


(108%); as a consequence, about 27% of ^o^Pb is missing with respect to (which has 

an absolute percentage solar abundance of 148%). Taking into consideration its still uncertain 
s-process contribution and the observational error in the determination of its solar abundance, we 
are missing about 10% of solar ^°^Pb at minimum. Note, however, that at odds with the s-only 
isotopes studied in this paper, this isotopes could receive a non negligible contribution from very 
low metallicity AGB stars (see Figure [B, which are not taken into account in our simplified GCE 
model. Thus, we can assume our ^'’^Pb production as a sort of lower limit. Concerning the test 
models previously discussed, we find that the absolute abundance of roughly scales as the 
*^°Sm one. On a relative scale, minor variations (<5%) are found in the majority of tests, apart 
from the Rotation case (-9%) and the Tail case (-1-22%). The latter could be a good candidate to 
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compensate the relative underproduction found in the Reference case. 

It is important to remark that, in addition to the uncertainties of AGB stellar models here discussed, 
other uncertainties may affect the predicted s-only distribution. As it is well known, AGB stars at 
various metallicity contribute d ifferently to the three s-process peaks. Thus, if the contribution 


from stars at large Z is favored (|Trippella et al. 


2014f) . a flatter relative s-only distribution may be 


found. Thus, the hypothesis on the existence of a LEPP process also relies on the uncertainties 
currently affecting Galactic chemical evolution models. We verified that an increase of the Star 
Formation Rate at all epochs leads to a faster increase of the ISM metallicity and, thus, to a larger 
contribution from metal-rich stars. As a consequence, we obtain a larger production of light s-only 
isotopes with respect to the heavy ones and, consequently, a flatter distribution. 

In conclusion, our full stellar evolutionary models coupled to a GCE model for the solar 
neighborhood does not necessarily require the need for a EEPP mechanism to be able to increase 
the Solar System s-only abundances in the range 96 < A < 124. However, owing to the 
uncertanties still affecting both stellar and Galactic chemical evolution model, we cannot a priori 
definitely rule out the existence of additional contributions to the Solar System s-only isotopes 
distribution. Note that the models presented in this paper cannot certify (or rule-out) the existence 
of a metal-poor primary EEPP, invoked to explain the abundances of a large group of light 
elements in low metallicity halo star which might be enriched by an r-process. 


We thank the referee for valuable comments and suggestions, that improved the quality 
of this paper. This work was supported by Italian Grants RBFR08549F-002 (FIRB 2008 
program), PRIN-MIUR 2012 ’’Nucleosynthesis in AGB stars: an integrated approach” project 
(20128PCN59) and from Spanish grants AYA2008-04211-C02-02 and AYA-2011-22460. 
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Table 1. Absolute percentage isotopic abunda nces with respect to the solar distribution (see text 


for details). Solar percentage errors (taken from 


Lodders et al. 


2009h are also reported (column 2). 


Isot. 

5o(%) 

Reference 

No IMS 

Rotation 

Tail 

Reimers 

Mloss AGB 

Nuclear 

SFR+10 

"OGe 

16 

17 

13 

17 

24 

17 

17 

15 

20 

^®Se 

7 

20 

16 

19 

30 

20 

21 

18 

24 

sokt 

20 

20 

16 

21 

32 

20 

20 

18 

25 

82Kr 

20 

26 

22 

25 

42 

27 

29 

26 

32 

86sr 

7 

52 

49 

47 

90 

53 

59 

52 

65 

«’Sr 

7 

51 

46 

45 

84 

52 

58 

52 

63 

‘^^Mo 

16 

no 

104 

68 

135 

113 

131 

104 

126 

IOOru 

6 

124 

118 

75 

151 

128 

148 

123 

143 

104pd 

11 

130 

124 

78 

157 

134 

155 

129 

148 

‘'°Cd 

7 

123 

117 

72 

144 

126 

146 

122 

140 

“'’Sn 

16 

112 

107 

65 

130 

115 

133 

109 

111 

‘22Te 

7 

126 

120 

70 

146 

130 

150 

117 

140 

i23Te 

7 

126 

120 

71 

149 

131 

149 

116 

141 

i24Te 

7 

135 

129 

74 

155 

140 

161 

124 

150 

‘28Xe 

20 

158 

152 

86 

181 

165 

190 

152 

111 

‘30Xe 

20 

154 

147 

82 

175 

160 

184 

148 

168 

‘34Ba 

18 

216 

208 

115 

257 

230 

255 

193 

236 

i36Ba 

18 

192 

185 

99 

217 

202 

230 

188 

207 

^^^Nd 

5 

170 

163 

80 

210 

184 

195 

174 

170 

‘^^Sm 

5 

159 

152 

75 

202 

173 

179 

152 

158 
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Table 1—Continued 


Isot. 

Se(%) 

Reference 

No IMS 

Rotation 

Tail 

Reimers 

Mloss AGE 

Nuclear 

SFR+10 

'^°Sm 

5 

148 

140 

68 

175 

159 

170 

145 

141 


14 

103 

97 

49 

115 

110 

117 

111 

102 

leoDy 

15 

139 

132 

65 

173 

151 

159 

134 

138 

ivoYb 

5 

154 

147 

72 

201 

167 

175 

150 

152 

‘^^Lu 

5 

183 

176 

83 

224 

186 

210 

189 

111 

176Hf 

5 

175 

174 

81 

225 

213 

199 

179 

111 

186os 

8 

111 

164 

81 

228 

189 

192 

170 

168 

192pt 

8 

128 

120 

62 

176 

140 

141 

124 

126 

198Hg 

20 

119 

112 

55 

149 

128 

134 

124 

116 

204pb 

7 

137 

130 

64 

160 

148 

152 

142 

133 







